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STOCHASTIC LANCHESTER-TYPE COMBAT MODELS I 



L. Billard 

Florida State University 



ABSTRACT 

There has been a lot of activity in recent years on 
the study of L an Chester -type combat models, especially from 
a deterministic standpoint. We consider some of these models 
in a stochastic framework and indicate how the appropriate 
deterministic model can be recast stochastically. Techniques 
for obtaining the corresponding solutions to the resultant 
differential-difference equations are discussed. These 
techniques are similar to those developed for use in other 
population process modelling situations such as epidemic 
theory and competition models. These are used to give the 
actual solution for the stochastic model of the original 
Lanchester model . 
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STOCHASTIC LANCHESTER-TYPE COMBAT MODELS I 



by 

L. Billard 

Florida State University 



1 . INTRODUCTION 

Since World War II there has been much activity on the 
study of combat models based on the original work of Lanchester 
(1914) . The so-called Lanchester models describe the situation 
in which two forces are in combat with each other with each side 
losing men or equipment (tanks) by attrition in accordance with 
some preassigned attrition law. These laws vary from the very 
simplest formulation used by Lanchester to quite complex cases. 

An extensive review of these different situation is presented 
in Taylor (1978) . 

Basically the model for the combat process can be 
described in terms of the sizes of the tv/o combat forces. Most 
of the work in the literature so far has confined its attention 
to studying the combat process in a deterministic framework. 

While such an approach can prove useful in providing broad guide- 
lines as to the behavior of a given combat situation, it is 
likely that a more accurate account can be obtained when the 
process is viewed stochastically. Therefore, in this paper it 
is shown how analogous stochastic combat models can be developed, 
and the appropriate solutions are derived. Specifically we 
show that a stochastic analogue of the Lanches ter-type combat 
models is nothing but a particular bivariate death process . 
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Furthermore if we wish to extend the simplest combat models 
to allow for reinforcements we will have a bivariate birth 
and death process . 

In the following section, we describe the general bivariate 
pure death process. This is applied to the particular case of 
the original Lanchester model in Section 3. In Section 4, we 

show how several other models can be formulated as a bivariate 
death process. These are models which generally have been 
considered in the literature from a deterministic approach. 

Some have been fully solved and others only partially solved 
(that is, for some finite subset of the governing parameters). 

We note that using the results of our Section 2, formal solu- 
tions can be found to the stochastic analogues in all cases. 

Then, in Section 5, we consider some models whose stochastic 
formulation is that of a bounded bivariate birth and death 
process. Finally we consider combat duration time in Section 6. 
Numerical results and comparisons for varying parameter values 
and battle force sizes will be presented in a companion paper. 
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2 . 



BIVARIATE DEATH MODEL 



Let B(t) and R(t) be the size of the Blue and Red 
forces at time t, respectively, and let particular state values 
be b and r, respectively. Typically, we may be interested 
in the number of men, or of tanks, or of both men and tanks, 
etc. We interpret "size" as the number of units (men, 
tanks, etc.) under discussion. Suppose that at time t = 0, 

B(0) = Bq and R(0) = Rq . 

The two forces are in combat with each other and units 
are lost according to various attrition laws, the particular 
attrition law depending on the actual modelling situation at 
hand. For general attrition rates a(b,r,t) and 3(b,r,t) for the 
Blue and Red forces, respectively, the general deterministic 
model is such that the following equations are satisfied: 



( 1 ) 



^ = -a(b,r,t) , 

^ = -6 (b,r, t) . 



When we view the process stochastically, that is, when 
B(t) and R(t) are taken to be random variables, the attrition 
terms a(b,r,t) and B(b,r,t) translate into infinitesimal 
transition probabilities according to the equations: 
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p{B(t+h) = b-1, R(t+h) = r|B(t) = b, R(t) = r} 

= a (b, r, t) h + o (h) , 

(2) P{B(t+h) = b, R(t+h) = r-l|B(t) = b, R(t) = r} 

= 3(b,r,t)h + o(h), 

P{two or more changes in (t,t+h)} = o(h) , 
and hence 

P{B(t+h) = b, R(t+h) = r|B(t) = b, R(t) = r) 

= 1 - {a(b,r,t) + 3(b,r,t)}h + o(h) , 

where o(h)/h = 0 . 

If we write 

Pb,r^^^ = P{B(t) = b, R(t) = r} , 

then the forward differential-difference equation governing this 
process is (there is a corresponding backward equation also) 

(3) ^ pj^ ^(t) = ~{a(b,r,t) + 3(b,r,t)} 

+ a(b+l,r,t) +3(b,r+l,t) Pt,,r+1^^^' 

for (b,r) € A = { (b, r) : 0 £ b £ Bq , 0 _< r £ Rq } , and where 
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Pb =0 whenever (b,r) / A. The initial conditions 

are 

‘ Pb,r<°> = ° 

As written here the transition probability generators 
a(b,r,t) and 3(b,r,t) are any (positive) function of the 
state of the process (b,r) and any function of time t. Fre- 
quently, the models of interest are such that these generators 
are time independent, that is, a(b,r,t) = a(b,r) and 
3(b,r,t) = 3(b,r). When this holds, the set of equations (3) 
is just the bivariate pure death process described by Billard 

and Kryscio (1977) . They then give the solution for p, (t) 

D , r 

for this equation for any generalized a(b,r) and 3(b,r) 
provided that 



(4) a(b^,rj^) + 3(b^,r^^) a(b2,r2) + 3(b2,r2) 



Thus , 



for (b^,r^) (b 2 ,r 2 ) ' 



(5) 




I I Cj^ (m,w|BQ ,R ) C 2 (m,w|b,r) exp{d (m,w) t} , 
m=b w=r 



where 



d(b,r) = -a(b,r) - 3(b,r) 

and where c^^ (m, w I B^ , ) and C 2 (ni,w|b,r) are given in Billard 

and Kryscio (1977, eqns . 10-11). 
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Most processes of practical interest are such that (4) 
holds naturally. In obtaining their solution, Billard and Kryscio 
(1977) essentially exploit and utilize a more general theorem 
given in Severe (1969a) . In the event that the restriction (4) 
does not hold, the solution can still be found by using the Severe 
theorem directly. Likewise, when the transition generators are 
functions of time t, direct use of Severe 's theorem yields the 
required solution to (2) . In any event, exploitation of the 
underlying structure which expresses itself as a partitioning 
scheme as used in Billard and Kryscio (1977) and Billard (1980) , 
has the advantageous effect of reducing considerably the degree 
of complexity that a first glance at Severe' s result suggests 
is involved. 
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3. LANCHESTER'S MODEL 



The simplest formulation for combat models is that of 
Lanchester (1914) where it is supposed that losses on each side 
are proportional to the size of the opposing force. This is the 
case of "aimed" fire. That is, 



a(b,r,t) = a(b,r) = 

and 

3(b,r,t) = S(b,r) = y 2 ^ • 



For convenience, we rescale so that a(b,r) = Ar and S(b,r) =b 
where now A is the relative effectiveness of the Red force to 
the Blue force. Then, the appropriate differential-difference 
equation governing this process is simply, from (3), 



(6) dt ^b r^^^ " -(Ar + b) 



b,r 



( t) + b 



^ Pb+l,r<P> 



for (b,r) € A and where we note that for boundary values 

(b,r) suitable adjustment is necessary. In matrix form, we may 
write (6) as 

= B(Pb,r<P>> 

where B is the matrix of coefficients (of transition generators) . 
Thus in the particular case that Bq = 4 , R^ = 2 and A = .8, 

B is given in Table 1 . 

The solution to (6) is given by (5) where it is readily 



verified that 
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(7a) Cj^(iti,w|Bq,Rq) 



B -m 

(Xw) ^ /(BQ-m)'. , 



w = Rq, 



T, . T • B»-m i, -irt 

R„-w+l i 2 j 0 0( ) 



0 

. I ■■■ . I 

\-w 



"0 ' 0 



^2"^1 "-Rq-w+I'^Rq-w 

• • • • • ^ 



<®0 - 1 - ^R(,-w> ' 



Rq-w 

X n 

V=1 



(B- + 1 - i ) 

0 V 

iBQ-m+l-i^+A (R q-w+1-v) } 



i -1 

V 



X n 



^v=^v-l 



{B_-m+l-j +A (R^-w+l-v) } 
0 -'v 0 



-1 



f w < Rq , 



with i. = 1 and ., = B„-m+l; and 

0 Rq~w+1 0 



(7b) C2(m,w|b,r) 



(-1)"' ^(Aw)"‘ V(in-b) ! 



r = w , 



‘^w-r+l ^2 i -i - m-b ^1 ^0, , /2"^1 

^ ••• I (-1) A (w-1) 



1 . 1 ~i 

w-r+1 w-r 



i ,=1 i =1 
w-1 1 



(A (w-r) + m - b} './{ A (w-r) +i -1} 

w~r 



w-r 

X n 

V=1 



i -1 

V 



(m-i +1) 

V 

(1-vA-i ) . 

V 3 =1 , 

V v-1 



n 



{a ( i-v) - 



-1 



r < w , 



with i = m-b+1. 

w-r+1 
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In matrix notation, the set of state probabilities may 



be written as 



^^b r ^ = C e (t) , 

where the vector e(t) has elements {exp(b^t)} with b^ 
being the ith diagonal element of B and the elements of the 
matrix C being determined from (7) . In our example, the 
corresponding C matrix is given in Table 2. Thus, for example, 

p{B(t) = 2, R(t) = 2} = 1.28e“^'^^ - 2.56e"'^*^^ + 1.283e"^'^^. 

Further details on how to calculate these quantities algorithmically 
will be discussed in the forthcoming paper dealing with numerical 
calculations and comparisons generally. 

We notice that the underlying structure for both c^^ ( * ) 
and C 2 (*) is very similar, and furthermore within each formula 
the structure for each i^ term is also similar. This allows 
computer calculation of these quantities to run without too much 
difficulty. This effect is especially apparent in the particular 
case that A = 1, that is, when the two forces are equally 
effective. In this case we find that 
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(m,w|BQ,RQ) • (Bq + Rq 



- m - w) 



B -m 
w , 



W = Rq , 



'Rq-w+1 



-Ro-w=l 



^2 1,-in 

J , ‘*0 



^2"^1 



Rq-w 

” ", ®o •" “- - ' 

V=1 



and 



C 2 (m,w|b, r) • (itH-w-b-r) 



, , , m-b m-b 
(-1) w , 



w 



^Rq-w+1 ^Rq-w 



w < Rq ; 



r = w , 



^w-r+1 ^2 . , 1 t~1a 1 1 

V V / ,xin+w-b-r 10 / 1,21 „ w-r-i-1 w-1 

I • * * z (“D w (w-1) r 

i =1 i =1 
w-r 1 



w-r 

X n ( m -I- 1 - i ) , 

r=l 



r < w 



Gye and Lewis (1976) have argued that it is quite reasonable that 
in the Battle of Trafalgar, A = 1. Other such situations can be 
easily visualized. 

Once expressions for Pj^ have been determined, other 

quantities of interest can be obtained. Suppose we are interested 
in the event that in time t the Blue force loses x units while 
the Red force loses none. Substitution in (5) and (7) yields 
the result 
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(8) P{B(t) = Bq - X, R(t) = Rq} 

= Pg _x,R ^ exp [- (Bq-x+RqA) t] (ARQ)^(l-e ^)^/x\ , 

X = 0 , ... , Bq . 

Likewise, if we are concerned with the probability that Red 
loses no units (regardless of the number of loses on the 
Blue side) up to time t, we have 

®0 

(9) P(R(t) = R„) = 

which can be determined from (8) . 

These results (8) and (9) could be likened to the model 
described by Gaver (1979) in which a force of initial size Bq 
attacks a bastion or stronghold of size R^ . Gaver assumed 

that the Red's stronghold is sufficient to guarantee no loss of 
Red units. He then considers the size of the Blue force at time t 
using always the deterministic approach. Gaver looks at two 
cases. One case assumes that the attrition rate of the Blue 
force is a result of unaimed fire by the Red force. Therefore, 
in the terminology of the present paper, the (stochastic) trans- 
ition rate is 

a(b,r,t) = p^(RQ/BQ)b, 

where is the attrition parameter. In the other case. Red 

emoloys aimed fire and hence the combat is modelled so that 

a (b, r, t) = p b 

3 . 
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for attrition parameter p^. In either case, Gaver's process 
is slightly different from the Lanchester results (8) and (9) 
since, although the quantity of interest includes the fact of no 
loses to the Red force, it is nevertheless possible probabilisti- 
cally that the Reds do lose some units by attrition. In Gaver's 
model it is assumed the Reds cannot lose any units whatsoever. 

Returning to the Lanchester model, we can easily obtain 
the expected number of units lost by time t by the Blue force 
given that the Red force has no loses. Let this random 
variable be denoted by X(t) . Then, 



B 



0 



Efx(t)} = ^ X P{B(t) = Bq - x|R(t) = Rq} 

x=0 



where 



“ I ^ Pr V r ^ Pb -X R 

x =0 x =0 ""'^0 



= AR„(l-e S(B--1)/S (B.) , 



'0 



0 



S (p) = \ {XR. (1-e ^)e^>Vy' • 



y=0 



If Bq is sufficiently large so that B(Bq) ~ S(Bq- 1), we have 



E{X(t) } ~ XR^d-e 



Hence, the eventual expected number of loses on the Blue side is 

lim E{X(t)} = XRq . 

t ->■ oo 

This result is what we would expect intuitively. 
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We can relate the stochastic and deterministic models 
through the expectations since the attrition rates are linear 
functions of the variables involved. Thus, if we multiply (6) 
throughout by b and sum over all (b,r) values, we have 

^ E{B(t) } = -AE{R(t) } . 

Likewise, multiplying by r and summing over (b,r) gives 

^ E{R(t) } = - E{B(t) } . 

Comparison with the deterministic equations shows that the 
expected values of the stochastic variables equals the solution 
of the deterministic equations. We note this correspondence 
does not necessarily hold when the attrition rates become non- 
linear in b and/or r. 
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4 . OTHER COMBAT MODELS 



In this section we briefly describe how other conibat models 
can be expressed stochastically. There have been many models 
proposed in the literature. An extensive and exhaustive review 
of such models has been made by Taylor (1978) to which we refer 
the reader for further details, elaboration, and justification. 

These models have by and large been studied from a deterministic 
viewpoint only. We confine outselves here to establishing certain 
stochastic analogues of some of these processes. Specifically, 
we obtain the appropriate format for the transition probability 
generators a(b,r,t) and 6(b,r,t). 

Once the generators have been established, it is then a 
relatively simple and direct procedure to solve the corresponding 
differential-difference equations as discussed earlier in Section 2. 
Obviously the simpler the form for a(b,r,t) and 3(b,r,t), the 
simpler the resulting solutions will be, such as we saw in the 
previous section for the original Lanchester model. Quite clearly 
the reduction achieved depends on the actual structure of these 
generators. When a(b,r,t) and 6(b,r,t) are linear functions 
of b or r, such reduction will be substantial and nice. How- 
ever, even the most complicated cases can still be solved since 
use of the Severo theorem is completely general without any con- 
fining restrictions as far as our models are concerned. 

Taylor and Brown (1976) establish what they call 
"variable-coefficient Lanchester- type equations of modern warfare" 



14 



where the attrition rate for each side is time-dependent . By 
analogy, we find the transition probability generators to be 
given by 

(10) a(b,r,t) = a ( t) r and 3(b,r,t) = b(t)b . 

If a ( t) = Yj^h(t) and b(t) = rescaling of time from 

t to h(t) produces generators 

a(b,r,t) = Yj_r and 6 (b,r,t) = Y 2 ^ ' 

which is the Lanchester model of Section 3. 

Range-dependent attrition rates were introduced by 
Bonder (1967) . Here, we have 

a(b,r,t) = -a(d)r and 6(b,r,t) = -6(d)b , 



where a(d) and B(d) are the attrition rates for the Blue and 
Red forces, respectively, dependent on the range or distance 
between the two forces. In some cases of tactical interest, 
a(d) and B(d) may be represented by the function 



and 



a(d) 



^1 

Pj_(l - d/d2) , 

0 , 



0 1 d _< d2 , 
d2 < d , 
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6(d) 



P2(l - d/d^) 



t 



0 _< d £ di 
dj^ < d , 

where d^^ and d 2 are the maximiun effective ranges of the Blue 
and Red forces, respectively, and where and ^2 

appropriate nonzero constants. If at t = 0, d £ M in (dj^,d 2 ) , 
that is, when combat begins both forces are within weaponary range 
of each other, and if there is no retreat, then this process re- 
duces to the Lanchester model of Section 3 where now 

= P^d - d/d3_^)^^, i = 1,2 . 

However, it is equally (perhaps more) reasonable to assume 
that the range parameter d is itself time dependent. Bonder 
(1967) takes 

d E d(t) = dp - vt 

where d^ is the initial range at t = 0 and v is the constant 

attack (closure) speed. Thus, the attrition rates a(d) and 

3 (d) are now time dependent. 

The range-dependent attrition rates of Bonder are a 

particular case of a more general power attrition rate class of 

models discussed by Taylor and Brown (1976) and Taylor and Comstock 

(1977) . This general class is such that the transition prob- 

^1 

ability generators assume the form a(b,r,t) = k^^(t + C) and 
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2 

3(b,r,t) = k 2 (t + C + A) , where A 0 is called the offset 
parameter and C 2. 0 is called the starting parameter. This 
accommodates the situation in which the forces have different 
maximum effective ranges and/or the battle begins within those 
maximum ranges. Taylor and Comstock (1977) give the solution 
to the deterministic model when there is no offset that is, 

A = 0, as well as for A > 0 when = 1 y^^ = 1, 

y 2 = 2. We note that the stochastic solution can be obtained 
for all y^^ and y 2 * 

A distinguishing characteristic of all the models considered 
so for is that the attrition rate of each force is proportional to 
the size of the other force and in no way depends on its own size. 

It is quite reasonable to expect that attrition could be propor- 
tional to the size of both opposing forces. This is especially 
so far unaimed fire (but can be equally argued for aimed fire) 
when it is realized that the probability that a unit on the Red 
force, say, will kill one unit (out of the remaining B(t) 
units) on tlie Blue force, increases (decreases) as the size B(t) 
increases (decreases) since there are more (fewer) units on which 
his fire may fall. In this context, the appropriate transition 
probability generators in a time independent situation would be 
a(b,r,t) = 3(b,r,t) = y 2 ^^ ’ note that here the 

generators are time independent so that the solution (5) holds 
directly. However, there is effectively a time dependency in 
that the actual size of each force is of course a function of time. 
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This model is mathematically equivalent to the Weiss (1963) 
predator prey model whose solution is given in Billard and 
Kryscio (1977) . 

Let us now consider the situation in which one or both 
sides has a supporting system which itself is not subject to 
attrition. Taylor and Parry (1975) considered a particular case 
of such a process which when presented as a stochastic process 
will have transition probability generators 

a(b,r,t) = -a(t)r - bj^(t)b , 

and 

3(b,r,t) = -b(t)b - a^^(t)r , 

where a(t) and b ( t) have the same interpretation as in (10) 
and where a^^(t) 0 and b^(t) 0 represent the attrition 

coefficients on the Red, and Blue, forces due to the supporting 
fire of the Blue, and Red, forces, respectively. Other situations 
come readily to mind such as a model in which the attrition from 
the supporting fire is also dependent on the size of its target 
(opposing force) . Then, we would have 

a(b,r,t) = -a(t)r - b^^(t)br , 

with 3(b,r,t) similarly defined if the same conditions prevail 
for the Red side. 
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Finally, Helinbold (1965) concerned himself with the 
situation in which grossly unequal force sizes are in combat. His 
modification of the Lanchester equations, when viewed stochastically, 
gives us 



and 



a(b,r,t) = -a(t) h(b/r)r 
3(b,r,t) = -b(t) h(r/b)b , 



where h(z) is interpreted as the effectiveness-modification 
factor. Helmbold imposed some conditions on the h(z). However, 
we note that for the stochastic model we are able to find the 
appropriate conditions without his restrictions 

Quite clearly the concept of supporting fire can be combined 
with Helmbold' s model. This has been done by Taylor (1976) in 
the particular case that h(z) = z for some constant c 0 . 
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5, BIRTH AND DEATH MODELS 

One feature of the models discussed to date is that only 
losses due to attrition occur to the combat forces, that is, they 
are pure death processes. However, there are some situations 
which may be adequately modelled as a birth and death process as, 
for example, when reinforcements are permitted. If the maximum 
number of possible reinforcements is known in advance, then we 
have a bounded birth and death process. 

As before, let B(t) and R(t) represent the size of 
the Blue and Red forces at time t, respectively; and let Bq = B(0) 
and Rq = R(0). Suppose that at time t there remain B*(t) and 
R*(t) Blue and Red reinforcements, respectively and suppose 
B*(0) = B* and R*(0) = R*, that is, B* and R* represent the 
maximum possible reinforcements available. Thus, after time t, 

B* - B*(t) actual reinforcements have been added to the Blue 
combat force, with R* - RgCt) similarly defined for the Red 
side. Let B* + B^ = B and R* + R^ = R. Let 

X(t) = (B*(t), B(t), R*(t), R(t) ) with realization x= (b*,b,r*,r) 
Then, for a completely general model, we have the following 
infinitesimal transition probabilities (corresponding to (2)) 

P{X(t+h) = (b*-l,b+l,r*,r) |x(t) = x} = Xj^(x,t)h + o(h) , 

P{X(t+h) = (b*,b,r*-l,r+l) |X(t) = x} = A 2 (X/t)h + o(h) , 

(11) P{X(t+h) = (b* ,b-l, r*, r) I X( t) = x} = Vij^(x,t)h + o(h) , 

P{X(t+h) = (b*,b,r*,r-l) lx(t) = x} = + o(h) , 

P{two or more changes in (t,t+h)} = o(h) , 
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and hence 



2 

P{X(t+h) = x|x(t) = x} = 1 - y {x.(x,t) + y.(x,t)}h + o(h) . 

i=l ^ ~ ^ 



If we write 



p (x, t) = P{X( t) = X} / 



the differential-difference equation governing the process is 
(compare (3)) 

^ 2 

(12) p(x,t) = - I {A^(x,t) + M^(x,t)} p(x,t) 

i=l 

+ (b*+l ,b-l , r*, r , t) p (b*+l ,b-l , r * , r , t) 

+ A2(b*,b,r*+l,r-l, t) p (b* , b, r*+l , r-1 , t) 

+ (b* ,b+l , r*, r) p (b* ,b+l , r*, r, t) 

+ ^ 2 (b*,b,r*,r+l) p (b* ,b , r* , r+1 , t) , 

for 



X € 8 



x:0 1 b* _< B*, 0 _< r* £ R* 
0<b <B, 0<r <R 



r 



and where p (x, t) E 0 whenever x / 8. The initial conditions 
are 



'® 0 '^ 0 ' ^0 



0 ) = 1 



p(x,0) = 0 for X ^ (B*,Bq,R*,Rq) . 
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The equation (12) can be solved by coit±)ining the techniques 
of Severo (1969a, 1969b) . When the transition generators are time 
independent, that is, (x, t) = A^(x) and vij^(x, t) = (x) , 

i = 1,2, an explicit solution for any arbitrary A^(x) and 
M^(x), i = 1/2, is given in Billard (1980). 

The vi£(x,t), i = 1,2, terms give the attrition rates 
and typically we would expect them to assume the same forms as 
the a(b,r,t) and 6(b,r,t), respectively, presented earlier. 

The terms Aj^(x, t), i = 1,2, represent the rates that reinforce- 
ments are added. Though written generally here, it is reasonable 
to expect that for models of practical interest, these generators 
will be a function of the size of the combat force itself. Thus, 
for example, the case where more reinforcements are brought in 
when the combat force itself becomes small could be reflected by 
a transition generator 



A^(x,t) = A(t)/b , 

or simply 

A^(x,t) = A/b . 

Recently, Gaver (1979) investigated the problem of infor- 
mation flow. In his model, individual units start off with un- 
aimed fire but then as information about the opposing force is 
gathered the fire becomes aimed fire. Thus, in the X(t) notation 
used above, B*(t) is the number of Blue units still using unaimed 
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fire at time t, B(t) is the number of Blue units using aimed 
fire at time t, and R*(t) and R(t) are the corresponding 
quantities for the Red force. If all units improve their fire 
to aimed fire before they themselves are killed, we have exactly 
the situation modelled in (11) and (12) . However, it is more 
realistic to assume that some of the units never learn from the 
information gathered, thus remaining with unaimed fire, and will 
themselves be lost by attrition. To model this situation, the 
infinitesimal transition probabilities (11) still apply in 
addition to 



p{X(t+h) = (b*-l ,b , r* , r ) |X(t) = x} = o(h) , 

and 

P{X(t+h) = (b*,b,r*-l,r) |x(t) = x} = o(h) , 



and where now 



P{X(t+h) = x|x(t) 



2 

X} = 1 - y {X. (x,t) + u. (x, t) + Y. (x,t) }h + o(h). 



Thus, for each i = 1,2, the (x, t) represent the 
information flow, the y^(x, t) represent the loss of aimed units 
from attrition, and the y 1) represent the loss of unaimed 
units from attrition. Suitable adjustment of equation (12) will 
give us the appropriate differential-difference equation for the 
process. This latter equation an then be solved using the method 
of Severo (1969a, 1969b) . 
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It is now possible to establish Gaver's (1979) model for 
information flow under mutual attrition in a stochastic frame- 
work. Thus, we have 

Aj^(x,t) = a^^b* and X 2 (X/t) = , 

where a and 3 represent the rate of conversion from an 
ua ua 

unaimed capacity to aimed capacity for the Blue and Red forces, 
respectively; 



yj^(x, t) = p^^r*b/B + 

and 

y 2 (Xft) = Y^^b*r/R + Y^^br/Cr + r*) , 

where represents the attrition rate of the unaimed Red 

forces against the aimed Blue forces, p^^ represents the attrition 
rate of the aimed Red forces against the aimed Blue forces, and 
Y and Y the corresponding attrition rates of the Blue 

forces against the Red forces; and 

Yj^(x/t) = p^^r*b*/B + p^^rb*/(b + b*) 

and 

Yo(x,t) = Y b*r*/R + y br*/(r + r*) , 

where p and p represent the rate of attrition to the unaimed 
^uu ^au 

Blue forces due to the unaimed and aimed fire, respectively, from 
the Red forces, and y Y similarly defined attrition 

UlU. 3\I 

rates of the Red forces. 
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6. COMBAT DURATION TIME 



The previous sections have been concerned primarily with 
techniques for obtaining the state probabilities of the underlying 
distribution. A knowledge of these probabilities allows most 
other quantities of interest to be derived. One such quantity 
is the combat duration time or the expected length of time for 
the battle. While this can certainly be obtained from our earlier 
results/ we present here a relatively easy straightforward but 
recursive method for its derivation. 

Suppose that at t = 0, the Blue force has Bq units 
and the red force has Rq units . Let us assume that the combat 
ends when the Blue force is reduced to a size of B' units and/or 
the Red force is reduced to a size of R' units. Hence, the 
permitted state space for (b,r) is now A' = {(b,r):B' £ £ Bq , 

R' £ £ Rq) instead of the A used in the previous section. 

(The adjustment from A to A ' in the results given earlier 
is trivial.) Let T(b,r) be the expected combat duration time 
from the time that there are b Blue units and r Red units 
remaining. Then, T(Bq,Rq) will be the overall combat duration 
time. For ease of illustration, let us suppose we have the simple 
Lanchester combat model where a(b,r,t) = Ar and S(b,r,t) = b. 
More general models are treated analogously. 

We first recall that the basic underlying process is 
Markovian. From (2), and basic properties of Markov processes, 
we see that when the system moves from the state (b,r) it moves 
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to the state (b-l,r) with probability Ar/ (b + Ar) or to 
the state (b,r-l) with probability b/(b + Ar) . Finally, before 
leaving the state (b,r) it will have stayed there for an average 
time of l/(b + Ar) . Therefore, we can write 

(13) T(b,r) = l/(b + Ar) + Ar/(b + Ar) T(b-l,r) + b/(b+Ar) T(b,r-1), 

for (b,r) € A'. In (13), T(b,r) =0 for (b,r) / A'. Clearly 

then, by starting at (b,r) = (B',R') and proceeding with 
(B',r), r = R'+l, ..., Rq, returning to (B'+l,r), r = R',...,Rq, 

and so on to (BQ,r), r = R Rq, we can find T(Bq,Rq). Note 
that in fact we have a matrix of T(b,r) values which give us 
the duration time from any intermediate stage (b,r) in addition 
to the overall duration T(Bq,Rq). 
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